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MARAAAAAAAAAAARALALLALALSLE SALES EELS LESSEE SERRE EEE EERE SEES EEE REESE TEER 


COPYRIGHT (c) 1978, 1980, 1982, 1984 B 
DIGITAL EQUIPMENT CORPORATION, MAYNARD. MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


® 
* 
® 
oY 
f 
THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION OF THE ABOVE CoPyni cut NOTICE. THIS SOFTWARE OR ANY OTHER = 
COPIES THEREOF MAY NOT GE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY © 
THER PERSON, NO TITLE TO AND OUNERSHIP OF THE SOFTWARE IS HEREBY * 
TRANSFERRED. : 
2 

THE INFORMATION IN THIS SOFTWARE 1S SUBJECT TO CHANGE WITHOUT NOTICE 
AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT * 
CORPORATION. . 
® 

® 

® 

te 

® 

® 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


® 
& 
® 
® 
& 
® 
® 
w 
-* 
-* 
® 
® 
:* 
-® 
-* 
-* 
-® 
-® 
-* 
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“FACILITY: General Utility Library 
ABSTRACT: 

Perform single precision floating point polynomial calculation. 
ENVIRONMENT: User Mode, AST Reentrant 


AUTHOR: Steven B. Lionel, CREATION DATE: 05-Oct-78 
MODIFIED BY: 
SBL, OF-Optores VERSION 1 

- Origina 


- Put HS number in standard format, with three digits in 
the edit number. JBS 16-NOV- 

Add "'_"’ to the PSECT ta ge ¥. Fy 
Minor code improvements. SBL 0 
Use local handler to insure tase “att exceptions other than those 
documented as being statuses are resignalled. SBL 25-Sept-1980 
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DECLARATIONS dba pate 11:09:46 tt 
9 g -SBTTL DECLARATIONS 
: ; INCLUDE FILES: 
§ SCHFDEF : 
f SSSDEF : 
09 23 ; EXTERNAL DECLARATIONS: 
00 61° .EXTRN LIB$SIG_TO_RET ; 
$0 $6 5 
4 ? 3 
00 4 : 
0000 65 
0000 66 ; 
0000 67 ; MACROS: 
0000 68 ; 
0000 4 
0000 0; 
0000 71 ; EQUATED SYMBOLS: 
0000 ie 3 
0000 7 
00000004 0000 74 arg = 4 : 
00000008 464 75 degree = § 3 
0000000C 000 76 coeff = 1 5 
00000010 0000 77 result = 16 ; 
0000 78 ; 
0000 79 ; 
0000 80 ; OWN STORAGE: 
0000 81; 
0000 82 
0000 83 ; 
44 Be 3; PSECT DECLARATIONS: 
00000000 86 -PSECT _LIBSCODE PIC, USR, CON, REL, LCL, 
9000 87 EXE, RD, NOWRT. LONG 


X/VMS Macro v04-00 
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Condition handling macros 
System symbols 


Library routine to convert 
a signal to an error return 
to caller of LIBSPOLYF. 

RO = signalled condition 


argument 
“ee of polynomial 
address of coefficient 


; result of polynomial 


table 


SHR, - 
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* -SBTTL LIBSPOLYF = Perform floating polynomial 
FUNCTIONAL DESCRIPTION: 


LIBSPOLYF provides the functionality of the VAX hardware 
instruction POLYF to high level language users. 


The third greress points to a table (array) of single 
precision floating point coefficients. The coefficient of 
the highest order term of the polynomial is pointed to 
by the table address operand, i.e. the first table element. 
The table is specified with lower order coefficients stored 
at increasing addresses. 
The evaluation is carried out by Horner's method, and the 
result is stored at the location pointed to by the fourth 
operand. The result computed is: 

if d = degree 

and x = arg 

result = CLOJ+x*(CC1J+x*(CC£2)+ ... x*C£dJ)) 


The unsigned word degree operand specifies the highest 
numbered coefficient to participate in the evaluation. 


For further detail, refer to the VAX-11/780 Architecture 
Handbook for the description of POLYF. 


CALLING SEQUENCE: 


status.wlc.v = LIBSPOLYF (arg.rf.r, degree.rw.r, coeff.rf.ra, 
result.wf.r) 


INPUT PARAMETERS: 
arg.rf.r - argument, ‘'x'’ in polynomial 
degree.rw.r degree of polynomial (GEQ 0) 
coeff.rf.ra - table of coefficients, floating 

IMPLICIT INPUTS: 

NONE 

OUTPUT PARAMETERS: 
result.wf.r - result of calculation 

IMPLICIT OUTPUTS: 
NONE 

FUNCTION VALUE: 
SS$_NORMAL 
SSS$_FLTOVF 


SS$~FLTUND 
55$~ROPRAND 


successful execution 

floating overflow 

te een under flow 7 
reserved operand, see VAX Architecture 
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8 Ve ; manual for more details 
0 13 ; SIDE EFFECTS: 
000 151: ALL other exceptions are signalled 
000 132; : : 
9000 184° 
400C¢ 009 155 -ENTRY LIBSPOLYF, “M<IV,R2,R3> ; Entry point, enable int. ovf. 
bo 138 3 and save R2, R3 
6D 15°AF 99E 809 138 MOVAB B“*HANDLER, (FP) ; Enable local handler to 
Babe 1) 3 exceptions 
Oc BC O08 BC 04 8C 55 0006 161 POLYF @arg(AP), - : perform polynomial 
000D 16¢ a@degree(AP), - 3; trap on exception to 
000D 16 acoef f (AP) 3; handler which will 
000D 164 3 unwind a return error 
Q00D 165 ; condition in RO to 
444 166 ; caller of LIBSPOLYF. 
10 Bc 50 £400 ones 198 MOVL RO, @result (AP) 3; return value 
50 01 9A oat 10 MOVZBL #1, RO 3 success status code 
04 0014 Ze RET : return 
0015 17 
0015 174 HANDLER: 
0000 0015 175 -WORD O 
ee 
0017 (178: If the exception is one of the documented exceptions for this routine, 
0017 179; call LIBSSIG_TO_RET to return it as a status. Otherwise, resignal. 
0017 180; Also, resignal Tf the depth is not zero. 
OO? dae 
50 O08 AC 00 0017 188 MOVL CHFSL_MCHARGLST(AP), RO ; Get mechanism vector address 
08 AO 05 OO1B 184 TSTL CHF SL_MCH_DEPTH(RO) : Is depth zero? 
41 12 OQO1E 185 BNEQ 90$ :; If not, resignal 
51 04 AC 00 0020 136 MOVL CHFSL_SIGARGLST(AP), R1 ; Get signal vector address 
50 04 A1 00 66 4 187 MOVL CHFSL_SIG_NAME(R1), RO ; Get signalled condition 
048C 8F 50 81 0028 133 CMPW =O, -#SS$_FLTOVF ; Compare conditions : 
gf 13 0020 189 BEQL 10$ ; If it matches, don’t resignal 
049C 8F 9 Bt ; ; 130 Haag ROE #SS$_FLTUND 
0454 8F 50 et B86 136 Haety ROE #SS$_ROPRAND 
04C4 BF 50 5 B23 133 is 8. #SS$_FLTUND_F 
04 Al 049C 8F * Ran 136 Nae #338_FLTUND. CHFSL_S!G_NAME(R1) ; Change fault code to trap code 
0484 8F 0 B} o31 138 5$: see OO #SS$_FLTOVF _F 
04 Al 048C 8F 3C 23 BY MOVZWL #SSS_FLTOVF, CHFS$L_SIG_NAME(R1) ; 
00000000'GF  6C FA 059 1 10$:  CALLG (AP)> G*LIB$SIG_TO"RET ; Return signal as a status 
50 0918 8F g $081 08 90$: MOVZWL #SS$_RESIGNAL, RO ; Resignal condition 
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Symbol table 66-SEP-1984 
ARG = 4 
CHF $L_MCHARGLST = : 
CHFSL-MCH_DEPTH s 
CHF SL_SIGARGLST s 4 
CHF 1G_N s 4 
OEF = 
DEGREE = 
HANDLER 15 R 0 
LISSPOLYF 0 00 RG 
LIBSSIG_TO_RET eeeerere 8 6X 
SULT = 010 
_FLTOVF = 48C 
SS$_FLTOVE _F = 4B4 
S$_FLTUND = 0000049C¢ 
SS_FLTUND _F = 000004C4 
SS$_RESIGNAL = 4's 4g bg 
SS$_ROPRAND = 00000454 
power ew econ es anaes + 
! Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
. ABS . 00000000 ¢ 0.) 00 ¢ 0.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
SABSS 00000000 ¢ 9.) 01 ¢ 1.) NOPIC USR CON ABS LCL NOSHR EXE RD WRT NOVEC BYTE 
_LIBSCODE 00000067 ¢ 103.) 02¢ 2.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
ee 
! Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 39 00:00:00.04 00:00:01.36 
Command processing 106 es SS 99: 89:06 -19 
Pass 189 88:80: 05°C 00:00:13.16 
Symbol table sort 0 00: B oR eS 00:00:01.18 
Pass 2 54 On°8 :00.5 00:00:01.3 
Symbol table output 4 00:00:06.0 00:00:00.0 
Psect synopsis output 2 i 00:00:00.0 
Cross-reference output 0 00:00:00.0 00:00:00.00 
Assembler run totals 386 00:00:04.0 00:00:19.16 


The working set Limit was 1050 pages. : 

21690 bytes (43 pages) of virtual memory were used to buffer the intermediate code. 

There were 30 pages of symbol table space allocated to hold 427 non-local and 3 local symbols. 
07 source Lines were read in Pass 1, produc ing 13 object records in Pass 2. 

mac 


pages of virtual memory were used to define ros. 
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' Macro library statistics a 


joensen wwe wmeeaancweaenr aan ames 


i 

Macro Library name Macros defined 
“S255S$DUAZB:(SYSLIBISTARLET.MLB;20 (iti(‘(‘;CS*~™S Te 
486 GETS were required to define 5 macros. 
There were no errors, warnings or information messages. 
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